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TRANSITION  TO  TURBULENCE  IN  CONSTANT-ACCELERATION  PIPE  FLOW 

INTRODUCTION 

Accelerating  shear  flows  and  other  time-dependent  laminar  and  turbulent 
shear  flows  are  encountered  in  many  aerod3mafflic  and  hydrod3mamic  applica¬ 
tions.  Examples  include  the  startup  of  a  closed-conduit  flow,  a  train  rapidly 
entering  a  tunnel,  flow  over  turbine  blades,  submarine  and  aircraft 
maneuverings ,  and  the  highly  transient  launch  of  a  torpedo. 

Even  in  simple  geometries  such  as  circular  pipes,  the  complex  nature  of 
unsteady  laminar,  turbulent,  and  transitional  flow  have  proven  to  be  very 
difficult  to  characterize  experimentally  or  solve  analytically.  Relatively 
little  theoretical  or  numerical  work  has  been  directed  at  this  problem. 

Only  recently  have  any  detailed  experiments  been  conducted  on  the  physics 
of  accelerating  flows.  The  most  comprehensive  study  was  a  previous  Naval 
Underwater  Systems  Center  (NUSC)  project,  which  was  conducted  over  a  3-year 
period  concluding  in  FY  1987  and  reported  in  reference  1.  This  was  the  first 
unsteady  flow  experiment  of  its  kind  conducted  at  NUSC. 

All  tests  of  that  project  were  conducted  in  the  NUSC  Unsteady  Flow  Loop 
Facility,  which  was  developed  and  built  coincident  with  the  project  and  funded 
mainly  with  NUSC  capital  asset  funds.  This  facility  is  capable  of  providing 
user-programmed  acceleration  of  the  mean  flow  over  a  large  range  of  velocities 
and  accelerations.  Details  of  this  facility  are  provided  in  reference  2. 

The  first  unsteady  flow  study  (reference  1)  used  a  5-cm-diameter  circular 
test  se-'tion  and  evaluated  the  effect  of  the  constant  acceleration  of  the  mean 
flow  on  many  of  the  prevalent  flow  parameters  including:  (1)  time-dependent 
velocity  profiles  and  turbulence  intensity  for  accelerations  starting  from 
rest  and  from  an  initially  turbulent  flow,  (2)  transient  wall  shear  stress, 
and  (3)  time  of  and  conditions  for  transition  from  laminar  to  turbulent  flow. 
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Earlier  experimental  studies  by  others  have  been  limited  in  i.he  types  of 
unsteady  flow  investigated.  In  general,  these  past  experiments  have 
investigated  either  pulsating  flow  (references  3,  4,  and  5),  a  suddenly 
applied  constant  pressure  gradient  (references  6  and  7),  or  a  suddenly  applied 
stepwise  change  in  flow  rate  (references  8,  9,  and  10).  In  short,  these 
earlier  experimental  studies  have  been  largely  limited  to  low  pipe  Reynolds 
numbers  and  in  the  flow  conditions  studied  and  the  types  of  data  acquired. 
Because  acceleration  and  the  time  rate  of  change  of  acceleration  changed 
continuously  over  any  one  run,  results  were  essentially  applicable  only  to  the 
specific  cases  studied. 

The  purpose  of  this  report  is  to  present  results  from  the  second  NUSC 
project  on  accelerating  flows,  which  was  conducted  during  FY  1988.  The  main 
objectives  of  this  second  project  were  to  investigate  the  details  of  the 
transition  process  and  to  validate  the  transition  correlation  parameters 
developed  under  the  first  project. 

The  details  of  the  transition  process  were  investigated  using  the 

existing  5-cm-diameter  test  section,  but  with  one  axial  location  supplemented 

with  additional  wall  shear  stress  sensors  in  an  effort  to  detect  the 

propagation  of  turbulence.  Validation  of  the  previously  developed  transition 

correlation  parameters  was  accomplished  by  conducting  constant  acceleration 

experiments  on  a  new  and  larger  9-cm-diameter  test  section.  Additionally, 

transition  for  mean  accelerations  lower  than  those  tested  previously  (i.e., 
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less  than  1.8  m/sec  )  was  investigated  on  the  9-cm-diameter  test  section. 


EXPERIMENTAL  FACILITY  AND  INSTRLT1ENTATI0N 


The  Unsteady  Flow  Loop  Facility  is  located  at  the  NUSC  Newport 
Laboratory.  The  design  operation  and  specifics  of  this  experimental  facility 
are  described  in  detail  in  reference  2.  A  general  overview  of  the  physical 
configuration  and  the  transient  flow  instrumentation  is  given  in  the  following 
paragraphs . 

The  layout  of  the  test  facility  is  shown  in  figure  1.  A  stainless  steel 
reservoir  supplies  a  150-hp  pump,  which  delivers  flow  to  the  test  section. 
Vibration  and  noise  isolators,  along  with  flow  straighteners  and  conditioners, 
are  provided  where  appropriate.  A  control  valve  installed  downstream  of  the 
test  section  is  part  of  the  feedback  control  system  that  provides  programmed 
acceleration  of  the  flow.  Reservoir  water  is  conditioned  to  maintain  its 
temperature  to  within  +1*C,  and  filters  remove  particles  as  small  as  0.5 
micron.  The  test  section  flow  rate  is  monitored  to  within  +1.0  percent  of 
reading  by  a  transient  electromagnetic  flowmeter,  whose  details  are  described 
in  reference  11.  This  fiowmeter  was  developed  as  part  of  the  flow  loop 
facility  design  project  since  there  are  no  other  known  commercially  available 
flowmeters  that  have  the  capability  of  making  accurate  transient  measurements. 

Figure  2  shows  the  5-cm-diameter  test  section  on  which  the  first  project 
was  conducted.  This  test  section  was  slightly  modified  for  the  first  phase  of 
the  present  tests.  The  test  section  consists  of  six  5-m-long  stainless  steel 
circular  sections  followed  by  a  30-cm  cast  acrylic  circular  section  through 
which  laser  Doppler  velocimeter  (LDV)  measurements  are  made.  There  is  a 
sensor  station  placed  at  the  end  of  the  six  stainless  steel  sections,  each  of 
which  consists  of  a  static  pressure  tap,  a  flush-mounted  wall  pressure  sensor, 
and  a  hot-film  wall  shear  stress  sensor.  These  sensor  stations  are  spaced  5-m 
apart . 

For  the  present  series  of  experiments,  the  acrylic  test  section  was 
fitted  with  four  hot-film  wall  shear  stress  sensors  positioned  near  the  LDV 
axial  measurement  station  as  shown  in  figure  2.  The  intent  of  the  instru¬ 
mented  acrylic  section  was  to  investigate  the  direction  of  propagation  of  the 
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Figure  1.  Flow  Loop  Facility  Layout 
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turbulence  at  transition  (i.e.,  does  the  turbulence  propagate  from  the  wall 
toward  the  pipe  centerline  or  vice-versa)  and  to  better  understand  the 
character  i'^*  ics  of  the  transitional  flow  field. 

ror  the  second  phase  of  the  present  project,  a  new  9-cra-dianieter  test 
section  was  built.  As  shown  in  figure  3,  this  circular  test  section  consists 
of  a  24.4-m-long  section  of  PVC  pipe  followed  by  a  6.0-m-long  cast  acrylic 
section.  The  cast  acrylic  section  was  used  mainly  for  LDV  measurements  while 
two  hot-film  wall  shear  stress  sensors  were  placed  in  the  PVC  section  at 
12.2-m  and  18.3-m  from  the  entrance  to  the  test  section  as  shown  in  figure  3. 
For  both  the  5-cm-diameter  and  9-cm-diameter  test  sections,  care  was  taken  to 
align  adjacent  pipe  sections  to  minimize  any  step  at  the  interface. 

For  the  9-cm-diameter  test  section,  the  control  valve  section  of  the  flow 
loop  facility  also  had  to  be  modified.  The  original  7.6-cm  Masoneilan  anti¬ 
cavitation,  anti-noise  cage  control  valve  was  replaced  with  a  larger  15.2-cm 
valve  to  accommodate  the  higher  flow  rates  of  the  larger  diameter  test 
section.  The  new  valve  was  a  Jamesbury  Corp.  high-performance  butterfly 
control  valve.  This  valve  provided  good  control  over  a  wide  range  of 
accelerations  and  velocities. 

For  all  test  runs  conducted  under  the  present  project,  the  control  system 

was  programmed  to  provide  constant  acceleration  from  rest,  leveling  off  to  a 

final  steady  mean  velocity,  which  was  sufficiently  above  the  transition 

velocity  so  it  would  not  affect  transition  itself.  Mean  flow  accelerations 
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ranged  from  0.2  to  11.15  m/sec  .  Transition  to  turbulence  was  monitored  by 
both  the  LDV,  focused  at  the  centerline  of  the  test  section,  and  the  wall 
shear  stress  sensors.  The  time  of  transition  was  obvious  in  the  abrupt  change 
in  the  value  of  the  measurements  at  transition. 

The  LDV  system  consisted  of  a  2-watt  argon-ion  laser  with  receiving 
optics  configured  in  a  forward-scattering  mode.  This  mode  was  used  to 
increase  the  valid  data  rate  from  the  counter  processor  to  a  value  between 
10,000  and  50,000  samples/sec.  Initial  tests  during  the  first  project  showed 
that  the  low  data  rates  typical  of  backscattering  techniques  (approximately 
400  samples/sec)  were  insufficient  to  follow  the  highly  transient 
instantaneous  velocities  of  these  experiments. 
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The  LDV  measuring  volume  was  0.93  mm  long  in  the  radial  direction  and 
0.161  mm  high.  Since  a  large  velocity  range  (more  than  30  to  1)  had  to  be 
accommodated  with  one  setting  of  the  LDV  counter  processor,  an  acousto-optic 
Bragg  cell  frequency  shifter  was  used  to  effectively  filter  out  the  LDV 
signal's  changing  pedestal. 

The  hot-film  wall  shear  stress  sensors  used  were  Thermal  Systems  Inc. 
(TSI)  1237-W  (flush  mounted).  The  hot-film  anemometer  used  with  the  wall 
shear  stress  sensors  was  a  TSI  model  #IFA  100.  All  data,  including  those  from 
the  transient  flowmeter,  wall  shear  stress  sensors,  and  the  LDV,  were 
collected  on  a  MASSCOMF  MC500  digital  data  acquisition  system  at  a  sampling 
rate  of  960  samples/sec.  A  specially  designed  interface  was  used  between  the 
MASSCOMP  system  and  the  LDV  counter  processor  to  permit  periodic  sampling  of 
the  LDV  data.  Details  of  the  instrumentation  and  the  test  procedures  are 
given  in  reference  1. 

It  should  be  noted  that  the  Apple  lie  computer,  which  was  part  of  the 
control  system  and  used  to  generate  the  transient  command  signal,  also  was 
used  to  trigger  the  start  of  the  data  acquisition  by  the  MASSCOMP.  The  Apple 
computer  then  waited  approximately  0.5  second  after  triggering  data 
acquisition  before  starting  the  acceleration  of  the  flow.  All  time  history 
curves  of  the  data  that  will  be  presented  in  this  report  will  display  this 
time  delay. 
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COMPARISON  WITH  THEORY 


The  transition  correlation  parameters  presented  for  the  9-cm-diameter 
test  section  data  require  the  calculation  of  boundary  layer  thickness  at  the 
time  of  transition.  As  shown  in  reference  1  for  the  transition  data  from  the 
5-cm-diameter  test  section,  the  velocity  profile  in  the  pipe  prior  to 
transition  agrees  with  the  profile  for  startup  of  laminar  pipe  flow  caused  by 
a  suddenly  applied  constant  pressure  gradient.  A  classic  analytic  solution 
for  that  case  was  given  in  reference  12.  The  solution  for  the  unsteady 
velocity  profile  U(r,t)  is  given  by: 


U(r.t) 

U 

max 


00 

=  (1  -  r*^)  -  I 


n=l 


8J-(X  r*) 
u  n _ 


K 

n  1  n 


t*) 


(1) 


2 

where  r*  =  r/R,  t*  =  ut/R  ,  and  X  are  the  zeros  of  the  Bessel 
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function  J^.  The  parameter  =  R  (-dP/dx)/4vi  is  the  final  maximum 

velocity  corresponding  to  steady  Poiseuille  pipe  flow.  For  t*  approaching 
unity,  U(r,t)  approaches  the  Poiseuille  flow;  for  t*  <  0.05,  the  centerline 
velocity  increases  linearly  with  time,  simulating  a  constant  acceleration 
startup. 


The  values  of  t*  at  transition  varied  from  0.000454  to  0.003654  for  the 
experiments  conducted  on  the  9-cm-diameter  test  section.  These  values  are 
small  enough  for  Szymanski's  solution  to  be  a  valid  approximation  for  the 
present  experiments.  Therefore,  equation  (1)  will  be  used  to  predict  the 
values  of  boundary  layer  thickness  at  transition  which,  in  turn,  will  be  used 
in  the  correlation  parameter-boundary  layer  thickness  Reynolds  number. 
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RESULTS 


In  this  section  of  the  report,  results  are  presented  first  for  the  data 
collected  from  the  5-cm-diameter  test  section  and  then  for  the  data  from  the 
9-cm-diameter  test  section.  The  9-cm  test  section  results  are  compared  with 
the  5-cm  test  section  results,  which  were  presented  in  reference  1. 

5-CM-DIAMETER  TEST  SECTION  RESULTS 

Prior  to  conducting  detailed  experiments  aimed  at  investigating  the 

physics  of  the  transition  process,  an  attempt  was  made  to  conduct  accelerating 

2 

flow  tests  from  rest  for  the  low  acceleration  range  (<1.81  m/sec  )  not 

tested  during  the  first  project.  Unfortunately,  the  control  valve  was 

inadequate  for  providing  the  very  fine  control  required  for  these 

accelerations;  a  smooth,  linear,  velocity-versus-time  curve  was  not  possible. 

However,  as  will  be  discussed  in  the  next  section,  the  control  system  for  the 

9-cm-diameter  test  section  was  able  to  provide  very  constant  accelerations  of 
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the  mean  flow  down  to  0.2  m/sec  . 

The  main  emphasis  during  testing  on  the  S-cm-diameter  test  section  was  to 
investigate,  in  detail,  the  physics  of  the  flow  field  during  the  transition 
process.  Specifically,  insight  into  the  generation  and  propagation  of 
turbulence,  as  observed  at  the  axial  location  where  LDV  measurements  were 
made,  was  desired. 

At  the  outset  of  the  present  project,  a  decision  was  made  to  limit  the 
investigation  to  whatever  information  could  be  obtained  by  using  only  the 
measurement  techniques  developed  and  used  during  the  previous  project. 
Therefore,  the  clear  acrylic  test  section  was  modified  to  include  four  hot-film 
wall  shear  stress  sensors  as  described  in  the  previous  section.  All  time 
histories  of  the  measurements  on  the  5-cm-diajneter  test  section  are  plotted 
with  960  points  per  second  of  data. 

A  typical  time  history  of  the  LDV  pipe  centerline  velocity  measurement 

2 

for  an  acceleration  of  1.86  m/sec  is  shown  in  figure  4.  The  abrupt  change 
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in  magnitude  followed  by  increased  fluctuations  about  the  mean  curve  are  a 
clear  indication  of  transition.  Figure  5  shows  the  corresponding  time  history 
for  shear  stress  sensors  6  through  10.  Outputs  for  sensors  1  through  5  are 
not  shown  since  they  are  similar  to  that  of  sensor  6,  which  was  at  the  end  of 
the  last  stainless  steel  section.  As  with  the  previous  project's  results 
reported  in  reference  1,  transition  occurred  at  approximately  the  same  time 
for  sensors  1  through  6  and  the  LDV. 


The  resulting  shear  sensor  data  for  stations  7  through  10  are  very 
unexpected.  As  shown,  these  four  stations  all  experienced  an  instability  that 
started  at  about  the  sane  time  (t  »  4  sec)  but  occurred  approximately  1.1 
seconds  prior  to  transition  to  turbulence  as  observed  at  the  LDV  and  the  other 
shear  stress  stations.  Station  10,  unlike  stations  7,  8,  and  9,  did  fully 
transition  abruptly  at  this  early  time  without  any  region  of  instability.  In 
figure  2,  it  is  shown  that  station  10  was  180*  from  station  8  while  stations 
7,  8,  and  9  were  placed  consecutively  along  the  flow  direction  and  spaced  4-cm 
apart.  Throughout  this  1.1-second  duration  of  instability,  the  shear  stress 
for  stations  7,  8,  and  9  experienced  high-level  fluctuations  that  were  larger 
than  the  levels  for  fully  turbulent  flow. 
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Figures  4  and  5  list  the  times  at  which  each  station  became  unstable  and 
then  fully  turbulent.  Station  10  became  unstable  (with  abrupt  transition) 
first.  At  0.037  second  later,  stations  7  and  8  both  became  unstable  at  the 
same  time  (within  the  0.001-second  resolution  of  the  data).  Station  9  became 
unstable  next,  0.004  second  later.  After  another  1.074  seconds,  a 
considerable  delay,  the  LDV  location  transitioned  to  turbulence.  This  is 
interesting  since  this  transition  at  the  pipe  centerline  was  1.115  seconds 
later  than  that  at  station  10,  which  is  at  the  same  axial  location.  Station  6 
then  transitioned  0.056  second  after  the  LDV.  Note  that  station  6  is  about 
25.6  cm  ahead  of  the  LDV.  Station  7  experienced  full  transition  0.033  second 
after  station  6  followed  by  station  8  after  another  0.014  second  and  then  by 
station  9  after  another  0.017  second. 

From  the  previous  measurements,  it  is  impossible  to  explain  the  flow 
phenomena  at  transition  except  to  say  that  there  is  very  little  indication  of 
propagation  effects  in  the  downstream  direction  and  there  appears  to  be  a 
substantial  amount  of  local  near-wall  effects.  The  only  propagation  type 
effect  can  be  seen  in  the  transition  at  stations  7,  8,  and  9.  Other  than 
that,  there  is  no  obvious  connection  between  the  various  stations. 

If  there  is  indeed  a  propagation  type  effect  from  stations  7  through  9, 
then  the  convection  velocity  of  the  transition  phenomenon  can  be  estimated 
from  the  time  information  to  be  about  2.5  m/sec.  Since  the  cross-sectional 
averaged  velocity  at  the  time  of  transition  for  these  locations  was  8.4  m/sec, 
the  convection  velocity  was  30  percent  of  the  cross-sectional  averaged 
velocity. 

The  extended  region  of  instability  was  never  observed  in  any  past 
experiments  when  the  wall  shear  stress  sensors  were  placed  only  in  the 
stainless  steel  sections.  It  is  interesting  to  note  that,  unlike  transition 
to  turbulence  at  stations  7,  8,  and  9,  the  time  of  initial  inst'bility  at 
these  stations  does  not  appear  to  have  any  axial  propagation  characteristics. 
This  tends  to  reduce  the  likelihood  that  the  reason  for  this  instability 
phenomenon  in  the  acrylic  section  is  that  the  joint  between  the  acrylic  and 
the  last  steel  section  acts  as  a  trip.  Also,  as  noted  previously,  great  care 
was  taken  to  minimize  any  steps  at  joints. 
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Another  fact  indicating  that  transition  in  accelerating  flows  is  not  a 
propagation  phenomenon  is  that  the  calculated  convection  velocity  of  transition 
from  station  6  to  7  is  inconsistent  with  that  between  stations  7  and  8  since 
it  is  7.8  m/sec  or  three  times  higher  than  that  between  stations  7  and  8. 

This  instability  phenomenon  occurred  consistently  for  repeat  test  runs  at 

the  same  acceleration  and  similarly  for  higher  accelerations  up  to  about  3 

2 

ra/sec  .  Beyond  that  acceleration,  the  region  of  instability  was  very  short, 
if  it  existed  at  all. 

Similar  data  to  those  just  presented,  but  for  an  acceleration  of  10.4 

2 

m/sec  ,  follow.  Figure  6  shows  the  time  history  of  the  LDV  output  for  this 
case,  while  figure  7  presents  the  shear  stress  sensor  output  for  stations  6 
through  10.  Again,  transition  occurred  at  approximately  the  same  time  for  the 
LDV  station  and  shear  stress  station  6.  Only  shear  stress  stations  8  and  9 
experienced  instability  prior  to  full  transition.  As  with  the  previous  test 
runs,  sensor  station  10  transitioned  earlier  than  all  the  other  locations  at  a 
time  of  0.67  second  before  the  LDV  and  a  substantial  0.18  second  prior  to 
station  9,  the  second  wall  location  to  become  unstable.  Station  8  then 
transitioned  0.014  second  after  station  9  first  experienced  some  Instability. 
Following  another  0.096  second,  station  7  transitioned.  The  LDV  and  station  6 
transitioned  0.384  and  0.394  second,  respectively,  after  station  7.  Obviously, 
this  "apparent"  transition  does  not  propagate  in  the  flow  direction.  In 
addition,  there  was  no  apparent  connection  between  stations  8,  10,  and  the 
LDV,  which  were  all  at  the  same  axial  location. 

Since  shear  sensors  7  through  10  were  fairly  distant  from  the  end  of  the 
acrylic  section,  it  was  difficult  to  ensure  that  the  sensors  were  flush  with 
the  wall.  A  protruding  sensor  could  locally  affect  the  apparent  transition. 
However,  all  the  sensors  tripped  relatively  early  to  turbulence,  reducing  the 
likelihood  of  protruding  sensors. 

The  previous  experiments  reported  in  reference  1  showed  that  the  velocity 
profile  do«fn  to  0.13  cm  from  the  wall  transitioned  virtually  instantaneously 
from  an  almost  uniform  profile  to  the  quasi-steady  one.  However,  from  the 
previously  mentioned  data  for  the  two  runs  presented,  it  appears  that  there  is 
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Figure  6.  vs  t;  X  »  10.4  m/sec  (5-cin  Test  Section) 

a  substantial  amount  of  near-wall  phenomena  occurring  locally  along  the  test 
section.  Since  the  boundary  layer  thickness  at  transition  ranged  from  &/R 
of  0.045  to  0.182,  it  is  possible  that,  at  some  axial  locations,  the  near-wall 
flow  acted  as  flat  plate  flow;  any  instabilities  remained  within  the  boundary 
layer  and  did  not  affect  the  core  flow.  Then  at  a  later  time,  full  transition 
occurred  virtually  globally  throughout  the  test  section  and  the  complete 
profile  became  turbulent. 

The  flow  field  in  the  vicinity  of  the  LDV  is  obviously  very  complex.  For 
a  complete  understanding  of  the  flow  field,  flow  visualization  is  required. 
Instantaneous  two-dimensional  flow  visualization  of  a  thin  slice  of  the  flow 
field  along  the  axial  direction  and  extending  from  stations  7  to  9  and  across 
the  diameter  would  be  ideal.  Potential  techniques  for  these  measurements  are 
now  being  developed  by  others  and  should  be  investigated  for  any  future  work 
in  this  area. 
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9-CM-DIAMETER  TEST  SECTION  RESULTS 

A  total  of  32  test  runs  were  conducted  on  the  9-cm-diaineter  test 
section.  Each  acceleration  run  was  started  from  rest  in  an  attempt  to 
determine  the  time  and  conditions  of  transition  to  turbulence.  All  figures  in 
this  section,  which  show  time  histories  of  the  measurements,  are  plotted  at 
120  points/second  of  data.  Note  that  in  these  time  histories  the  time 
t  =  0.0  second  is  at  the  start  of  data  acquisition  and  that  the  actual 
acceleration  begins  approximately  0.5  second  later. 

Figures  8  through  10  show  the  resulting  linear  cross-sectional  averaged 

velocity  versus  time  for  three  different  test  runs  at  accelerations  of  0.2, 

2 

6.3,  and  11.15  m/sec  ,  respectively.  These  runs  cover  the  range  of 

2 

accelerations  tested,  which  varied  from  0.2  to  11.15  m/sec  .  The  high 
linearity  of  these  curves  verifies  that  constant  acceleration  was  indeed 
obtained.  On  each  of  the  figures,  the  time  at  which  transition  occurred  (as 
evidence  by  the  wall  shear  stress  sensors  and  the  LDV)  is  noted.  It  is 
obvious  that  transition  to  turbulence  did  not  affect  the  acceleration. 

It  was  fortunate  that  the  flow  loop  facility,  as  configured  for  the  9-cm- 
diameter  test  section,  could  provide  substantially  lower  accelerations  than 
was  possible  with  the  5-cm  section.  As  shown  in  the  following  discussion  of 
results,  the  new  data  add  substantially  to  the  data  base  from  the  5-cm  tests. 

Figures  11  through  13  are  the  time  histories  of  the  output  from  the  two 

wall  shear  stress  sensors  for  the  three  runs  of  figures  8  through  10  (i.e, 

2 

accelerations  of  0.2,  6.3,  and  11.15  m/sec  ,  respectively).  Since  only  the 
time  of  transition  was  of  interest  and  not  the  absolute  magnitude  of  the  wall 
shear  stress,  the  sensors  were  not  calibrated.  The  time  of  transition  is 
clear  where  the  sensor  output  increases  almost  instantaneously. 

For  each  of  the  32  runs,  the  time  of  transition  was  essentially  global, 
that  is,  transition  occurred  abruptly  and  at  virtually  the  same  instance  in 
time  as  indicated  by  the  various  measurements.  Transition  times  generally 
differed  by  less  than  3  percent  with  a  few  runs  at  a  maximum  of  10  percent 
when  based  on  the  time  from  the  start  of  the  acceleration. 
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Table  I  lists  the  correlation  parameters  calculated  for  the  32  runs. 

Test  runs  are  presented  in  order  of  increasing  acceleration  and  include 

=  final  steady  state  velocity,  T  =  temperature,  v  =  kinematic 

viscosity,  U  ^  =  cross-sectional  mean  velocity  at  transition,  t  = 

m.tr  ^  ’  eff,tr 

time  of  transition  from  the  start  of  the  transient  (not  the  start  of  data 

acquisition),  and  =  nondimens ional  time  at  transition  based  on 

t  .  Note  that  t  ^  is  the  effective  time  of  transition  as 

ei t , tr  et  t , tr 

calculated  by  t  ,,  =  U  /X  to  account  for  the  delay  between  the  start 

6  L  X  f  L  IT  tn  f  w  F 

of  data  acquisition  and  the  transient. 


As  in  reference  1,  the  transition  correlation  parameters  are  Re^  = 

0  y  t  IT 

pipe  Reynolds  number,  K  =  local  acceleration  parameter  at  transition, 

f  C  IT 

and  RCf  ^  =  boundary  layer  thickness  Reynolds  number  at  transition,  where 

o ,  tr 


Re„  ^  =  D  *  U  /v, 

D,tr  m,tr 


(2) 


Re.  „  =  S,  * 

6,tr  tr 


m,  tr 


(3) 


K 

a 


(v/U 


3 

ro,  tr 


)  *  (dU  /dt). 

OD 


(4) 


The  first  correlation  of  the  transition  data  is  shown  in  figure  14  where 
acceleration  is  plotted  against  Re-  .  The  curve  for  the  9-cm  data  is 

U  f  wIT 

shown  along  with  that  for  the  5-cm  data.  The  data  for  each  individual  test 
section  do  follow  a  trend  of  increasing  Re-  with  increasing  acceleration, 

1/  f  c  ^ 

but  the  data  for  each  do  not  agree  with  each  other.  These  parameters, 

therefore,  are  obviously  inadequate  as  a  transition  correlation,  but  the  data 

show  that  all  of  the  values  of  Re„  ^  are  extremely  high  (up  to  Re-  ^  of 
^  D,tr  j  o  r  D,tr 

1.11  X  10°)  relative  to  the  steady-state  transition  value  of  approximately 
2000  for  pipe  flow. 


Figure  15  shows  Re_  ^  vs  t*^  .  It  can  be  seen  that  neither  of  the 
®  D,tr  tr 

two  parameters  are  constant  at  the  transition  point  of  an  accelerating  flow. 
However,  as  stated  in  reference  1  for  the  5-cm-diameter  test  section,  these 
two  parameters  do  seem  to  correlate  with  each  other.  At  the  lower  Re^ 
values,  the  correlation  between  the  5-cm  and  9-cm  data  is  reasonably  good;  the 
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Table  1.  Transition  Parameters;  9-cm-Diameter  Test  Section 
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Figure  14.  Re^j  vs  X  at  Transition 
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9-cin  section  data  being  only  slightly  lower  but  following  the  same  trend  as 
the  5-cm  data.  For  the  higher  Re^^  tr^^®D  tr  ^  ^  *  iO^),  the  values  of 
Re_  approach  a  constant  value  of  approximately  1.1  x  10^.  Again,  even 
though  the  parameters  do  seem  to  correlate  well  between  test  sections, 

Re^  does  not  appear  to  be  an  adequate  correlation  parameter  since  it  does 
not  include  realistic  size  or  length  scale  information.  That  is,  as  shown  in 
reference  1,  the  boundary  layer  at  the  wall  at  the  time  of  transition  is  very 
small  and  does  not  extend  to  the  pipe  centerline,  thereby  making  the  pipe 
Reynolds  number  questionable. 


In  reference  1,  an  attempt  was  made  to  develop  a  correlation  parameter 
for  local  acceleration,  which  was  the  analog  of  the  convective  acceleration 
parameter,  K  =  (v/U  )(dU/dx),  which  is  commonly  used  for  relaminarization. 

The  local  acceleration  parameter  that  was  developed  was  K  ,  as  given  by 

di 

equation  (3).  A  detailed  discussion  of  the  parameter  K  is  given  in  reference 

1. 


The  curve  of  K  vs  Re_  is  plotted  in  figure  16.  For  the  lower 

Si  U  f  C  IT 

Re^  where  the  9-cm  and  5-cm  data  overlap,  good  agreement  for  the  two 
cases  is  evident.  For  the  larger  Re^j  the  9-cm  data  follow  the  same 
trend  as  observed  for  the  5-cm  data  in  that  the  value  of  K  asymptotically 

3i 

approaches  a  lower  value  as  Re_,  increases.  For  the  9-cm  data,  the  value 

-8  _g 

is  approximately  0.7  x  10  .  The  data  for  K  <  U.O  x  10  are  replotted 

d 

in  figure  17  to  give  a  better  indication  of  the  trend.  It  is  obvious  that 
this  curve  appears  to  be  a  reasonable  choice  for  correlating  transition  data; 
however,  it  does  not  appear  that  the  lower  value  of  K  has  as  yet  been 

Si 

reached . 


The  last  attempt  described  in  reference  1  to  arrive  at  a  reasonable 

correlation  was  to  plot  Re-  .  versus  acceleration.  This  correlation 

6  9 1  r 

includes  realistic  length  scale  information  in  Re^  Figure  18  is  a 

plot  of  these  variables  for  both  the  previous  5-cm  data  and  the  present  9~cm 

2 

data.  For  the  range  of  accelerations,  1.8  to  11.8  m/sec  (where  data  were 
obtained  on  both  test  sections),  excellent  correlation  between  the  two  sets  of 
data  are  evident. 
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For  the  acceleration  range  above  3  m/sec  ,  the  value  of  Re^  is, 

o ,  t  r 

as  desired,  approximately  constant.  The  mean  value  of  Re^  for  t.^e  5-cm 

o ,  tr 

data  is  2^984  +  7.8  percent,  while  for  the  9-cm  data  the  value  is  26511  +  9.3 

2  ~ 
percent.  For  the  50  test  runs  above  3  m/sec  from  both  the  5-cm  and  9-cm 


test  section  data,  the  mean  value  of  Re^,  ^  is  25,595  z  10.4  percent. 

6  9  tr 

The  scatter  about  the  mean  values  were  calculated  at  the  95-percent  (two 
standard  deviations)  confidence  level. 


2 

For  the  lower  accelerations,  0.2  to  1.8  m/sec  ,  obtained  on  the 

9-cm-diameter  test  section,  the  value  of  Re^  decreases  with  decreasing 

acceleration.  In  fact,  the  shape  of  the  mean  curve  for  this  range  is  very 

similar  to  that  predicted  through  an  extrapolation  in  reference  3.  This  trend 

with  decreasing  acceleration  is  seen  to  asymptotically  approach  the  steady 

state  (zero  acceleration)  transition  Re^  value  of  approximately  2000  for 

flat  plate  flow.  It  is  interesting  to  note  that  for  the  lowest  acceleration 
2 

tested,  0.2  m/sec  ,  which  is  a  mere  1/49  of  the  acceleration  due  to  gravity, 
the  value  of  Re^  is  approximately  10,500,  still  five  times  higher  than  the 
zero  acceleration  value. 


This  trend  of  decreasing  Re^  is  substantiated  by  the  single  experiment 

of  van  de  Sande  et  al.  (reference  7).  They  report  transition  during  the 

acceleration  phase  of  startup  flow  in  a  pipe.  In  their  experiment,  they  used 

a  quick-opening  valve  and  a  constant-head  tank  to  generate  accelerating  flow 

in  a  5-cm-diameter  water  pipe.  Their  experiments  were  limited  to  very  low 

accelerations.  For  one  of  their  runs,  transition  occurred  at  t  -t-  4.2  sec, 

corresponding  to  a  mean  velocity  of  1.15  m/sec  and  a  Reynolds  number 

Re_  ^  +  57500.  From  figure  8  in  reference  7,  one  can  estimate  that  the 

D,tr  *  2 

very  low  acceleration  at  the  transition  point  was  about  0.1  m/sec^.  The 

estimated  value  of  Re^  a  8600  is  given  in  figure  18  and  is  shown  to  follow 

the  trend  of  the  present  9-cm  data. 

An  interesting  observation  was  made  for  three  of  the  test  runs  conducted 

on  the  9-cffl-diaineter  test  section.  As  shown  in  figures  19  and  20  for  the 

2 

accelerations  of  1.25  and  7.1  m/sec  ,  respectively,  the  flow  at  wall  shear 
stress  sensor  number  2  tripped  to  turbulence  earlier  than  expected  but  then 
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relaminarized  prior  to  the  final  transition  at  the  expected  time  coincident 

with  transition  at  sensor  number  1,  In  figure  21,  for  an  acceleration  of  5.65 
2 

m/sec  ,  a  similar  phenomenon  occurred  but,  this  time,  sensor  number  1 
exhibited  the  premature  transition. 

Some  local  instability  is  the  likely  cause  for  the  initial  transition. 

However,  from  the  limited  wall  shear  stress  measurements,  it  can  only  be 

assumed  that  the  reversion  back  to  a  laminar  type  signal  is  due  to  the  passing 

of  some  sort  of  turbulent  puff,  which  was  confined  to  a  relatively  short 

length  of  the  test  section.  By  using  the  known  mean  velocity  at  the  beginning 

and  end  of  the  first  turbulent  region,  the  puff's  length  was  estimated  to  be 

2 

1.3  m  for  the  1.25  m/sec  acceleration  of  figure  19,  2.0  m  for  the  7.10 

2  2 
m/sec  acceleration  of  figure  20,  and  1.2  m  for  the  5.65  m/sec 

acceleration  of  figure  21. 

It  is  unlikely  that  relaminarization  due  to  the  stabilizing  effect  of  the 
acceleration  occurs.  Some  justification  for  that  statement  can  be  obtained 
from  previous  experiments  reported  in  reference  1.  There,  relaminarization 
was  not  observed  for  an  initially  low  velocity  turbulent  flow  undergoing 
constant  acceleration.  This  interesting  phenomenon  would  be  a  good  topic  for 
further  investigation. 

It  should  be  noted  that  for  the  three  runs  that  exhibited  the  two 
transition  points,  the  t*'s  presented  earlier  were  calculated  using  the  time 
of  the  final  transition  point. 
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CONCLUSIONS 


A  unique  series  of  experiments  was  conducted  on  the  laminar-to-turbulent 
transition  in  constant-acceleration  pipe  flow.  The  major  conclusions 
resulting  from  this  study  are: 

1.  The  main  effect  of  mean  acceleration  on  the  flow  field  is  to 
stabilize  the  flow  to  such  a  great  extent  that  transition  to  turbulence  is 
considerably  delayed. 

2.  Pipe  Reynolds  numbers  at  the  point  of  transition  in  accelerating 
flows  have  been  observed  to  be  as  high  l.l  x  10^,  substantially  higher  than 

3 

the  steady-state,  zero-acceleration  value  of  approximately  2  x  10  . 

3.  The  correlation  parameters  for  predicting  transition  in 
constant-acceleration  pipe  flow,  namely,  K  and  Re;,  have  been  validated 
as  being  reasonable  and  accurate. 

4.  The  details  of  the  flow  field  at  the  time  of  transition  are  very 
complex.  It  appears  that  some  local  effects  lead  to  local  instabilities  and 
possibly  transition  that  are  confined  to  the  near-wall  region. 

5.  A  comprehensive  understanding  of  the  flow  field  is  impossible  without 
a  global  investigative  technique  such  as  laser-speckle  velocimetry  flow 
visualization. 
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